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Relativistic effects are investigated in nuclear matter calculations employing renormalized low- 
momentum nucleon-nucleon (NN) interactions. It is demonstrated that the relativistic effects cure a 
problem of non-relativistic low-momentum interactions, which fail to reproduce saturation of nuclear 
matter. Including relativistic effects, one already obtains saturation in a Hartree-Fock calculation. 
Brueckner-Hartree-Fock calculations lead to a further improvement of the saturation properties. 
The results are rather insensitive to the realistic NN interaction on which they are based. 
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A fundamental question in nuclear physics is the role 
played by relativistic effects. At first sight, only rather 
moderate relativistic effects are expected, since the ve- 
locity of nucleons in nuclear matter reaches less than one 
fourth of the light velocity and the kinetic and poten- 
tial energies are small as compared to the rest mass of 
the nucleons. However, considering a meson-exchange 
model for the NN interaction the weak single-particle po- 
tential occurs as a sum of a very attractive scalar field 
and a repulsive vector field, each of them of a size as 
large as several hundred MeV. This leads to a modifi- 
cation of the Dirac spinors for the nucleons in the nu- 
clear medium and can be expressed in terms of an effec- 
tive Dirac mass, which tends to decrease with increasing 
density [H, Si. The matrix elements of the NN inter- 
action in the nuclear medium should be evaluated us- 
ing these dressed Dirac spinors and therefore depend on 
the density of the matter considered. These relativis- 
tic effects lead to a successful microscopic description 
of the saturation properties of nuclear matter in Dirac- 
Brueckner-Hartree-Fock (DBHF) approach^, 0, @| with- 
out a need to include many-nucleon forces, which are 
required in non-relativistic investigations @. 

Apart from the division between relativistic and non- 
relativistic approaches, one can distinguish between in- 
vestigations which are based on phenomenological inter- 
actions and those which are based on realistic NN in- 
teractions. The parameters of the former interactions 
have been adjusted to describe properties of isospin sym- 
metric nuclear matter and of nuclei in the valley of 
f] stability. These phenomenological models, such as 
Gogny 0, Skyrme [g], or relativistic mean field (RMF) 
approaches provide a simple description of the mean 
field in terms of local single-particle densities. However, a 
disadvantage is that their predictive power may be rather 
limited, in particular for highly isospin asymmetric nu- 
clear matter and nuclei far away from the line of f3 sta- 
bility. The study of exactly these nuclei is of high inter- 
est with the forthcoming new generation of radioactive 
beam facilities like the future GSI facility FAIR in Ger- 
many and the RIA facility planned in the United States 
of America. 



One the other hand, the microscopic approaches 
based on realistic NN interactions, like the Bonn 
interactions have a high predictive power, which 
should give one confidence when the model is used in 
extreme cases like in a highly isospin asymmetric nu- 
clear environment. These high-precision free-space NN 
interactions are adjusted to describe the experimental 
data of the NN interaction. However, the strong short 
range and tensor components of such realistic interac- 
tions make it inevitable to employ non-perturbative ap- 
proximation schemes for the solution of the many-body 
problem. Therefore, calculations with such interactions 
are restricted to very light nuclei due to the dramatic 
increase of configurations with an increasing number of 
nucleons. Furthermore, such sophisticated calculations 
will not become feasible for heavier nuclei or the nuclear 
structures in the neutron star crust in the near future. 

A possible way out of this problem is to restrict the nu- 
clear structure calculation to the low momentum compo- 
nents by separating the low momentum and high momen- 
tum components of a realistic NN interaction by means 
of renormalization techniques 0, El, [i2- If the cutoff 
A is appropriately chosen, i.e. around A = 2 fm^^, the 
resulting low momentum interaction Viowk will still de- 
scribe the experimental data of the NN interaction up 
to the pion threshold. Moreover, a very attractive aspect 
is that this Viowk interaction turns out to be independent 
of the underlying realistic interaction V. 

Using Viowk iu a calculation of nuclear matter, one ob- 
tains a binding energy per nucleon increasing with den- 
sity in a monotonic way H, 13, 14\, unless three-body 
forces are added Thus, the emergence of a satu- 

ration point is prevented in symmetric nuclear matter. 
Attempts have been made to cure this problem and re- 
cently Siu et al.[l3] suggested to increase the cutoff, ac- 
count for the contribution of pphh ring diagrams and a 
modification of the NN interaction, the so-called Brown- 
Rho scaling [16j]. Increasing the cutoff parameter leads to 
the necessity to account for correlation effects {pphh ring 
diagrams) in a non-perturbative way and thereby one of 
the nice features of the low momentum interaction, that 
they can be treated in a perturbative way, is lost. 
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Therefore, in the present work, we want to follow a 
different route: We include the relativistic effects dis- 
cussed above and evaluate for each density the matrix 
elements of the bare NN interactions using the Dirac 
spinors, which are appropriate for this density. In this 
work, the medium properties of the nucleons, which are 
used to dress these Dirac spinors, are obtained from the 
EoS presented in Refs. [5|]. From this density-dependent 
bare interaction we can than deduce a renormalized Viowk 
interaction, which will depend on the density, using the 
standard techniques. 

For the construction of a low momentum potential, 
it is necessary to separate the low momentum and high 
momentum components of realistic interactions. In our 
work, the unitary-model-operator approach (UMOA) 
17| is used to disentangle these parts. As in all well- 
known model space techniques, we define the projection 
operators P, which projects onto the low momentum sub- 
space, and Q, which projects onto the complement of 
this subspace, the high momentum subspace. Further- 
more, these operators P and Q satisfy, of course, the 
usual relations hke P + Q = 1, ~ P, Q"^ — Q, and 
PQ = = QP. The aim of the unitary-model-operator 
approach is now to define a unitary transformation U in 
such a way that the transformed Hamiltonian does not 
couple P and Q, which means 



QU-^HUP = 



(1) 



has to be fulfilled. It leads to an effective Hamiltonian 



eff, 



(2) 



which contains a starting Hamiltonian ho describing the 
one-body part of the two-body system and an effective 
interaction Veff. This effective interaction is defined in 
terms of this unitary transformation as 



Veff = U-\ho + V)U-h„, 



(3) 



with V representing the bare NN interaction. The uni- 
tary operator U can be expressed as 



J7 = (1 + - cjt)(l + ujuj'' + Lu^fuj)-^^^ 



(4) 



with an operator to fulfilling the relation uj = QujP such 
that — = 0. This operator uj can be obtained by 
first solving the two-body eigenvalue equation 

{ho + V)\<i>k)^Ek\^k) 



(5) 



and afterwards defining the matrix elements of to us- 
ing the eigenstates |$p) having the largest overlap with 
the low momentum space, the P-space Uj. Next, the 
effective interaction Veff is calculated as described in 



111 Il8l |. In this way one obtains the effective Hamil- 



tonian of Eq. (21) , which contains the effective interaction 



Veff. The eigenvalues, which are obtained by diagonal- 
izing this effective Hamiltonian in the P-space, are iden- 
tical to those, which are obtained in the diagonalization 
of the original Hamiltonian H = ho + V in the complete 
space. 

Next, this model-space scheme can be applied to the 
effective two-nucleon problem by considering for the basis 
states of the two-nucleon system the states identified by 
the relative momentum, its modulus and the correspond- 
ing partial wave. For a given partial wave, the states with 
a relative momentum smaller than a cutoff A are identi- 
fied as the states of the P-space. Therefore, applying the 
technique described above leads to the normal effective 
interaction Viowk. 

In order to address the question which role correlations 
beyond the simple mean-field or HF approximation play 
for the density dependent Viowk {p) interaction, the effects 
of NN correlations are taken into account by means of 
the BHF approximation. In the BHF approach, one re- 
places the bare NN interaction V, or in our case the 
effective interaction Vef / , used in the HF approximation 
by the G-matrix which obeys the Bethe-Goldstone equa- 
tion. In the two-particle center of mass frame, it takes 
for an effective interaction Vef / with cutoff A the form 

.A 

G{k',k,ek)eff = Veff{k',k)+ / q^dqVeffik',q) 

Jo 

7. ^f—-G{q,k,eu), (6) 

2efc - 2eq + irj 

where with i = k,q are single-particle energies. Fur- 
thermore, the Pauli operator Qp prevents scattering to 
occupied states and, therefore, restricts the intermediate 
states to particle states with momenta q which are above 
the Fermi energy. In Eq. ([SI), we also take into account 
that Vef f is designed for a model space with relative mo- 
menta smaller than A. Therefore, the integral in Eq. ( [6|) 
is restricted to momenta q below the cutoff A. In con- 
trast, there is no upper integration limit in case of bare 
NN interactions. 

Results for the energy per nucleon of symmetric nu- 
clear matter as a function of the density p obtained from 
HF and BHF calculations are displayed in Fig. [TJ In all 
our calculations we use a cutoff A of 2 fm~^ and the Bonn 
aQI interaction has been employed for the calculations 
displayed in this figure. The conventional method, which 
ignores the medium modifications of the Dirac spinors, 
leads to results with the well known features: The HF 
calculation using Viowk does not exhibit a minimum in 
the energy as a function of the density. This absence of 
saturation is one of the major problems of calculations 
for nuclear matter employing Viowk- This problem can- 
not be cured by the inclusion of correlations beyond the 
HF approximation, e.g. by means of the BHF approxi- 
mation. 

These rather weak effects of the NN correlations be- 



3 



> 




0.1 0.2 0.3 0.4 

-3 

Density p [fm ] 



FIG. 1: Binding energy per nucleon of symmetric nuclear 
matter of a HF calculation (dotted line) and of a BHF calcu- 
lation (dashed line) employing a density dependent Viowk(p) 
interaction. Furthermore, comparison is made with a HF and 
a BHF calculation using a standard Vio-wk interaction, i.e. a 
Viowk calculated at density p = 0. . . 



mean-field calculations using bare NN interactions p. 
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FIG. 2: Binding energy of nuclear matter as a function 
of density from HF calculations employing various underly- 
ing potentials for the density dependent effective interaction 

Vlowk{p)- 



yond mean-field can be explained by the fact that only 
intermediate two-particle states with momenta below the 
cutoff A can be taken into account in Eq. ([H]), which is not 
the case for bare NN interactions. The NN correlation 
effects even vanish at high densities due to the lack of 
phase space for these correlations and energies resulting 
from BHF calculations approach those determined in the 
HF approximation. 

The inclusion of relativistic effects by calculating the 
underlying Bonn A interaction in terms of dressed Dirac 
spinors, however, drastically changes this picture. The 
HF calculation for the density dependent Viowk{p) in- 
teraction already leads to a minimum in the energy as a 
function of the density. Furthermore, this HF calculation 
yields a saturation density in the neighborhood of the em- 
pirical region of saturation, but it yields too little bind- 
ing energy, i.e. the energy is about -11 MeV per nucleon. 
This result for the binding energy can be improved by 
taking into account the effects of NN correlations within 
the BHF approximation. Hence, the BHF calculation for 
the density dependent Viowk[p) interaction yields about 
2 MeV more binding in the region of saturation than the 
corresponding HF calculation due to the additional corre- 
lations from NN states. In short, the mean-field calcula- 
tions using a density dependent Viowk[p) can already lead 
to reasonable results and the NN correlations beyond 
mean-field are rather weak, in particular compared to 
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FIG. 3: Binding energy of nuclear matter as a function of 
density from BHF calculations employing various underly- 
ing potentials for the density dependent effective interaction 

Vlowkip)- 
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An attractive feature of the use of a standard Viowk 
interaction in HF or BHF calculations is that the cor- 
responding results are rather insensitive to the realis- 
tic interaction model on which the Viowk interaction is 
based. Therefore, we want to explore to which extent this 
insensitivity applies in nuclear matter calculations em- 
ploying the density dependent Viowk{p) interaction. The 
HF calculations for density dependent Viowkip) interac- 
tions based on different realistic interactions qualitatively 
shows the same behavior, although quantitatively some 
difference in the binding energy strength exists as can be 
seen in Fig. [2] The HF calculation using Bonn A as un- 
derlying interaction yields about 2.2 MeV stronger bind- 
ing in the saturation point than the one using Bonn C. 
The same observations and conclusions about the sensi- 
tivity on the underlying realistic interaction can be made 
for the BHF calculations in Fig. [S] This dependence of 
the calculated energy on the underlying NN interaction 
is much weaker than the corresponding difference of 8 
MeV per nucleon obtained in non-relativistic BHF calcu- 
lations using directly Bonn A and C potentials 0]. The 
model-dependence of the relativistic Viowk interactions, 
however, seems to be slightly stronger than the corre- 
sponding model-dependence in non-relativistic calcula- 
tions using the modern potentials, which fit the NN data 
with high precision. Here, one must keep in mind, how- 
ever, that the potentials Bonn A, B and C did not fit the 
phase-shifts that accurately and indeed a similar model- 
dependence is also obtained applying the renormalisation 
scheme to these potentials in the conventional way. 

In summary, we study relativistic effects of nuclear 
matter calculations employing an effective low momen- 
tum interaction. These effective low momentum interac- 
tions are very interesting, since calculations with realistic 
bare NN interactions are restricted to very light nuclei 
and not possible for heavier nuclei or for the neutron 
star crust in the near future. However, one of the ma- 
jor problems of these calculations employing an effective 
low momentum interaction is the absence of saturation, 
unless three-body forces are added. Therefore, relativis- 
tic effects are introduced by dressing the Dirac spinors 
of the underlying realistic interaction on which the low 
momentum interaction is based. In this way a density 
dependent effective interaction Viowkip) is obtained. Em- 
ploying this density dependent Viowkip) interaction, one 
already obtains a saturation point in the HF calculation. 



The effects of the NN correlations introduced by the 
BHF approximation, although rather weak, leads to a 
further improvement of the saturation properties. Hence, 
this leads to the nice feature that HF calculations using a 
density dependent Viowkip) can lead to reasonable results 
and NN correlations beyond mean-field are rather weak. 
This opens the door for studying finite nuclei in calcu- 
lations which are based on a realistic NN interaction, 
treating correlation effects in a perturbative manner. 
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